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Abstract 

Using the Matsubara formalism, we consider the massive {Xip^)o vector A^-component model in 
the large N limit, the system being confined between two infinite paralell planes. We investigate 
the behaviour of the coupling constant as a function of the separation L between the planes. For 
the Wick-ordered model in = 3 we are able to give an exact formula to the L-dependence 
of the coupling constant. For the non- Wick-ordered model we indicate how expressions for the 
coupling constant and the mass can be obtained for arbitrary dimension D in the small-L regime. 



Closed exact formulas for the L-dependent renormalized coupling constant and mass are obtained 
in D = 3 and their behaviours as functions of L are displayed. We are also able to obtain in 
generic dimension D, an equation for the critical value of L corresponding to a second order phase 
transition in terms of the Riemann 2;eia-f unction. In D = 3 a renormalization is done and an 
explicit formula for the critical L is given. 
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1 Introduction 



Translationally non-invariant field theories are of great interest in several branches of theoretical 
physics, particularly in the study of confined systems, or systems undergoing changes of state in- 
duced by space inhomogeneities. For instance, in the Ginzburg-Landau theory of phase transitions 
applied to the study of type II superconductors, the external magnetic field, constant in some fixed 
spacial direction breaks translational invariance introducing important changes in the Feynman 
rules [jl]], 0, 10. Another important situation is found in the finite temperature field theoretical 
formalism of statistical mechanics, in particular the Matsubara imaginary time formalism, which 
introduces a breaking of translational invariance along the imaginary time (temperature) axis, 
leading to corresponding modifications in the Feynman rules. In this context, in recent papers 
(01) [0]) 0) the behaviour of coupling constants on the temperature has been investigated in con- 
nection to stability and phase transitions. In the finite temperature behaviour of the X(f% model 
was investigated at one-loop approximation. One of the results is that the thermal squared mass 
(in the disordered phase) is a positive increasing function of the temperature, while the thermal 
correction to the renormalized coupling constant is negative and increasing in modulus with the 
temperature. This could be interpreted as the possibility of the vanishing of the thermal coupling 
constant at some temperature and its change of sign afterwards. Of course we can not conclude 
from this one-loop result that there is a breaking of stability at some temperature, since we can 
not be sure that this peculiar behaviour will be preserved when higher order loop's contributions 
are considered (see below). 
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Questions concerning stability and the existence of phase transitions may also be raised if we 
consider the behaviour of field theories as function of spacial boundaries. In this situation, for 
Euclidean field theories, time (temperature) and the spacial coordinates are on exactly the same 
footing and the Matsubara formalism applies as well for the breaking of invariance along any 
one of the spacial directions. The possibility of vanishing of the coupling constant as mentioned 
above for temperature, or more generally the existence of phase transitions, would be in this case 
associated to values of some spacial parameters describing the breaking of translational invariance, 
for instance the distance L between planes confining the system. In this case, in a very similar way 
as it has been remarked in the context of finite temperature field theory in ||^, two-loop corrections 
in the small L regime are positive and tend to compensate the one-loop lowering of the value of 
the L-dependent renormalized coupling constant. Similar questions could be raised in general 
about the existence of phase transitions induced by spacial boundaries separating different states 
of the system. One of the most striking examples of this kind of situation is the spacial separation 
between the Abrikosov lattice of vortices and the other phases in type II superconductors [0. In 
another related domain of investigation, there are systems that present defects, created for instance 
in the process of crystal growth by some prepared circumstances (domain walls). At the level of 
effective field theories, in many cases this can be modelized by considering a Dirac fermionic field 
whose mass changes of sign as it crosses the defect, what means that the domain wall plays the 
role of a critical boundary, separating two different states of the system (see for instance 
and other references therein). In any case the question of the behaviour of physical parameters 
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in field theories (masses or coupling constants) on spacial boundaries or inhomogeneities, is still 
largely open and deserves further investigation. 

In order to shed some light on the problem, in this article we examine a model that allows a 
non-perturbative approach. We consider the Euclidean vector A'"-component X(fi% theory at leading 
order in the system being confined between two infinite parallel planes a distance L apart from 
one another. In this case we are able to obtain for the Wick-ordered model an exact formula for 
the renormalized coupling constant behaviour as a function of the given spacial boundaries. For 
the non- Wick-ordered model we obtain for the coupling constant and the mass, formulas valid in 
the small-L regime for arbitrary dimension D. In dimension D = 3 we obtain closed formulas 
for the renormalized L-dependent coupling constant and mass and plots of these parameters as 
functions of L are given. We are also able to obtain in generic dimension D, an equation for 
the critical value of L corresponding to a second order phase transition in terms of the Riemann 
zeta- function. In D = 3 a renormalization is done and an explicit formula for the critical L is 
given. 

2 The Wick-ordered model 

We consider the model described by the Lagrangian density. 




(1) 
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in Euclidean D-dimensional space, where Aq is the bare couphng constant, m is the physical mass 
(see below) and summation over repeated "colour" indices a is assumed. We note that Wick- 
ordering of products of fields makes unnecessary to perform an explicit mass renormalization at 
the order at the one-loop approximation. The tadpoles are suppressed by Wick-ordering. Thus 
we consider in the following only the zero external momenta four-point function, which we take as 
the basic object for our definition of the renormalized coupling constant. To simplify the notation 
in the following we drop out the colour indices, summation over them being understood in field 
products. 

The four-point function at leading order in is given by the sum of all the diagrams of the 
type depicted in Fig.l. We consider the system confined between two parallel planes a distance 
L apart from one another, normal to the x -axis. We use Cartesian coordinates r = [x, z), where 
z is a (£) — l)-dimensional vector, with corresponding momenta k = (k^, q), q being a (D — 1)- 
dimensional vector in momenta space. Then since the x-dependence of the field is confined to a 
segment of length L, we are allowed to perform the Matsubara replacements, 

fdk, 1+2? 2n7r 

Performing the sum over all diagrams indicated in Fig.l we get for the L -dependent renormalized 
four-point function at zero external momenta, 

where T^{D,L) corresponds to the single bubble diagram in Fig.l, after subtraction of the ultra- 



violet divergent parts. 

To obtain an expression for Ti{D,L), we briefly recall some one-loop results described in ||^ 
adapted to our present situation. These results have been obtained by the concurrent use of dimen- 
sional and zeta-function analytic regularizations. In order to deal with dimensionless quantities 
in the regularization procedures, we introduce parameters = m^/47r^/i^, {L^Y = a^^, g = 
(A/Stt^), (v^o//^) = (po, where ipo is the normalized vacuum expectation value of the field (the 
classical field) and /i is a mass scale. In terms of these parameters the boundary- dependent {L- 
dependent) one-loop contribution to the effective potential can be written in the form, 

oo /'_lNs +00 „ ^D-lh. 



or, using a well-known dimensional regularization formula ||T0[ , 



L) = /i^v^^ /(D, s)g^<Pl^zf{s - — — ; a), (5) 

s=l 

where f{D, s) is a function proportional to r(s — ■^^) and Zf{s — a) is one of the Epstein- 
Hurwitz zeta -functions defined by, 

2 '^'^ 

Z^(u;ai,...,ax) = E {ainl + ... + axrij^ + c^y, (6) 

ni,...,nK=—oo 

valid for Re{u) > K/2 (in our case Re{s) > D/2). The Epstein-Hurwitz 2;eta-f unction can be 
extended to the whole complex s-plane and we obtain, after some manipulations the one-loop 
correction to the effective potential. 



-| oo 

U,{D,L) = -Y.g'<Pl'h{D,s) 

^ s=l 



-n / J-y.,IIL N J^-^a ^ — III, . IJ . , 



(7) 



where h{d, S) = 2d/2-s-\^2s-d/2 ^r{s+l) ^"^^ Bessel functions of the third kind. 

Coming back to the situation under study, we see that the single bubble function is 
the coefficient of the fourth power of the field (s = 2) in the above equation. We can write S(Z}, L) 
in the form, 

j:{d,l) = h{d)-g{d,l), (8) 



where 



and 



I 2S-D/2 D-4(_iY+l jj 
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oo r 



2 (2vr)^/2 



-1 D/2-2 



KD_2{mnL). (10) 



.nL_ 

^From the properties of Bessel functions it can be seen from Eq.(p!OD that for any dimension D 
G{D, L) satisfies the conditions, 

limG(D,L)=0 lim G(L>, L) ^ oo. (11) 

We also conclude from Bessel functions properties that G(-D, L) is always positive for any values 
of D and L. 

Let us define the renormalized coupling constant A in the absence of boundaries as, 

1a= limrg)(0,L). (12) 

^From Eqs.(0) and (H) we get, 

^^'^ = ^ 1-Ao [H{D)-G{D,L)] ^^^^ 



and from (0), (0) and (|T3|), 



Combining (|T^) and ([TsD , we obtain the L -dependent renormalized coupling constant, 

rg)(0.L)4A(AL) = i^^^^. (15) 

The above procedure corresponds on perturbative grounds, to sum up all the chains of bubble 
graphs in Fig.l with all possible combinations of and counterterms -^S\ = at the vertices. 
It is just the resummation of all the perturbative contributions including the counterterms from 
the chains of bubbles and the subtraction of the divergent (polar) parts, written in compact form. 
These subtractions are to be performed even in the case of odd dimension, where no poles of 
F-functions are present. 



An exact result can be obtained in dimension D = 3. From the relationships ||12|| , 



^n+i W = K.{z) = y^e-^ (16) 

we obtain from Eqs.([T3|) and (PUD, after summing a geometric series, the exact expression for the 
coupling constant in the large N limit, 

3. L) = /""TV (17) 

^ ^ 8m7r(e™-^ - 1) + 3A ^ ^ 

where the subscript W is used to indicate explicitly Wick-ordering. A plot of Xw{D = 3,L) is 
given in Fig. 2. 
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3 The Non- Wick-ordered model 

The effect of suppression of Wick-ordering is tliat tlie renormalized mass can not be taken as tlie 
coefficient of tlie term ipa'^a in tlie Lagrangian, as we have done in the preceding section. We must 
take a L -corrected renormalized mass which satisfies the gap equation in the ordered phase, 



m%(L) = H ; — rrrTrr 



Ln 



Kn_,{mR{L)nL), (18) 



where m stands for the constant mass parameter in the Lagrangian density (|1|), which corresponds 
to the renormalized mass in the absence of boundaries (free space) for the Wick-ordered model. 
To obtain the L -dependent coupling constant the constant mass parameter m should be replaced 
in Eqs(|TB|) and ( [TI]| ) by the L-corrected mass mij(L) and the resulting system of equations should 
be solved with respect to L. Under these conditions exact closed expressions are almost impossible 
to obtain, since it would require a procedure equivalent to solving exactly the Dyson-Schwinger 
equations. 

Nevertheless results can be obtained, both analytic and numerical. In this sense in the next 
we indicate how to generalize to our case in arbitrary space dimension D, some results already 
obtained in refs. and fl^ for D = 4 in the context of thermal field theory. We take an integral 



representation for the Bessel function in Eq. ( p^Sf ) |T2 



M^) = ^^(|)'^ / ^"^*(^' - ^r^^'dt, (19) 
valid for Re{v) > — | and \arg{z)\ < |. Using this representation in Eq.(|^), after some straight- 



forward calculations the L-dependent renormalized mass can be written in the form 



D-S 



raid) = rn' + '-^±^F(D),nr(L) f -^rrr (^TT^? " 1 I (20) 



where F{D) = ^ g„i p/a-ix - When D is odd, the power is an integer and Newton binomial 



theorem gives a direct way for evaluating the integral in Eg. (pO]) , implying in an algebraic equation 
for m\{L). When D is even the expansion of — l) ^ yields an infinite power series. 

We obtain, 



^l^L) =m' + -^^j^ E 9{D, k){mn{L)Lr ^ (21) 

^Viv JmB.{L)L e — i 

where 5'(-D, fc) = l)*^C o-3 and the C's are the generalized binomial coefficients for arbitrary 
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power. For small values of k the integral that appears in the above equation is the well known 
Debye integral 



oo 



A-,-)=/ ^-^ = Ee--a:"(i + ^ + ... + ^), (22) 



which is valid for a; > and n > 1. For k > which corresponds to n < 1 in the preceding 
equation, the exponent of r in eq.(|^) becomes negative and the integral is undefined. In this case 
for small values of L, a generalization to negative odd powers of the argument of the integrand in 
the Debye integral can be done using the methods of ref. |T3| and the integral has the expansion 
{u = mFi{L)L and n = D — 3 — 2k), 

J{u,n)=l dr— - = - V — p -Bn\nu + an (23) 

Ju e--l ^^^^^ q\ q-n n\ 
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where Bk are the Bernoulh numbers and a„ is a constant. Replacing the above equations in 
Eq.(21) we have in the small L regime the following expression, 



ml{L) =m' + ^Y^Z? [^(^' D) + B{L, D)] , (24) 



where, 



and 



A{L,D)= 9{D,k){mR{L)Lf''I{mn{L)L,D-3-2k) (25) 

A;=0 



fc=oo 

B{L,D)= J2 9{D,k)mR{L)LY^J{mR{L)L,D-3-2k) (26) 

The equations above give a non-perturbative expression for the L-corrected renormalized mass in 
the small L regime, for even dimensional Euclidean space. Closed expressions are very difficult 
to obtain in arbitrary dimension, but for L sufficiently small, the series in Eq . (|26|) and the repre- 
sentations (??), (p3|) for the functions /(m/j(L)L, D — 3 — 2k) and J{mfi{L)L, D — 3 — 2k) could 
be truncated, giving an approximate algebraic equation for the L-corrected mass. The resulting 



expression should be injected in eqs.(|T0D and ( p!5D to give the L-corrected couphng constant in the 



small L regime. We will not perform these manipulations here. Instead, we will in the following 
obtain some exact results in D = 3. 

In dimension D = 3, using Eqs.(|l6D the sum over n in Eq.(|l^) can be performed exactly. We 
obtain a closed transcendent equation satisfied by the L-corrected renormalized mass, 

KiL) =m'- log(l - e-«(^)^), (27) 
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or, in the large N limit, 



ml{L) = - A log(l - e-'"^^^)^). (28) 



The large N renormalized mass is plotted as a function of L in Fig.3 

Also, as in the preceding section, we have in = 3 an exact expression for the coupling 
constant in the large limit as a function of the renormalized L-dependent mass, 

Xn{D - 3,L) - 8^^^(^)(^„,(^), _ i) + (29) 
In Fig. 4 we plot in dimension D = 3 the L-corrected coupling constant for the non- Wick-ordered 
model using jointly Eg s . (1281) and (|29|) in the same scale we have used to plot the Wick-ordered 
coupling constant. We see that the behaviour of the coupling constant is quite different from 
the monotonically increasing in behaviour of the renormalized squared mass and also from the 
behaviour of the coupling constant in the Wick-ordered case. The non- Wick-ordered coupling 
constant slightly decreases for decreasing values of L until some minimum value and then starts 
to increase. In the Wick-ordered model the coupling constant tends monotonically to zero as L 
goes to zero, while in the non- Wick ordered model it has a non vanishing value even for very small 
values of L. In fact, numerical analysis of the solution of Eq. (^Sj) shows that mji{L)L and 
m\{L)L — > cxD as L — > and, therefore, the L-corrected non- Wick-ordered coupling constant has 
a non- vanishing value at L = 0, which is equal to the free space value A. As a general conclusion it 
can be said that for the non- Wick-ordered model the L-dependent renormalized coupling constant 
slightly departs to lower values, from the free space coupling constant. Furthermore this departure 
is smaller for smaller values of A (in units of m). 
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In what concerns the mass behaviour, let us refer to Eq.(|T8|) and remember that the Bessel 
function K is positive for real positive values of the argument and decreases for increasing values 
of the argument. For space dimension D > 2 the correction in L to the squared mass is positive 
and the L-dependent squared mass is a monotonically increasing function of If we start in 
the disordered phase with a negative squared mass, the model exhibits spontaneous symmetry 
breaking of the 0{N) symmetry to 0{N — 1), but for a sufficiently small critical value of L the 
symmetry is restored. The critical value of L, is defined as the value of L for which the inverse 
squared correlation length , y^o), vanishes in the gap equation. 



where (po is the normalized vacuum expectation value of the field (different from zero in the ordered 
phase). In the neighbourhood of the critical point ipo vanishes and the gap equation reduces to 
Eq.([l8|). If we limit ourselves to the small L regime, we may use an asymptotic formula for small 
values of the argument of the Bessel functions in Eq. (p^. Then taking m/j(L) = in the resulting 
equation it is not difficult to obtain the critical value of L in Euclidean space dimension D {D > 2), 



. ,31) 



where 



9{D) = -^T{^-l)aD-2), (32) 

4:71- 2 

({D — 2) being the Riemann zeta-function. For D = 3 the zeta -function in g{D) has a pole and 
a renormalization procedure implies that the quantity L^'"^ is proportional to the regular part of 
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the analytic continuation of the 2;eta-function in the neighbourhood of the pole. It is not difficult 
to show that in this case the critical value of L is given by the solution of the equation, 

where the parameter \x corresponds to a fixed renormalization scale. The result above generalizes 



to a phase transition associated to a spacial boundary, estimates and numerical simulations in |T6 
and W%. 



4 Concluding Remarks 

Taking the Wick-ordering, which eliminates all contributions from the tadpoles, we decouple in 
some sense the boundary behaviour of the coupling constant from the mass boundary behaviour. 
Wick-ordering is a useful and simplifying procedure in the applications of field theory to par- 
ticle physics, but the same is not necessarily true in applications of field theory to investigate 
critical phenomena, where the contribution from tadpoles could be physically significant. As a 
consequence of the suppression of Wick-ordering the boundary behaviour of the coupling con- 
stant is sensibly modified with respect to the monotonic behaviour in the Wick-ordered case (see 
comments following Eq.(^)). As a general qualitative conclusion it can be said that, for the 
non- Wick-ordered model in the perturbative regime in A, the L-dependent renormalized coupling 
constant only slightly departs from the free space value A and even so, for small values of L. 
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Figure captions 

Fig.l - Typical diagram contributing to the four-point function at leading order in To each 
vertex there is a factor and for each single bubble a colour circulation factor N 

Fig. 2 - Renormalized coupling constant (in units of j^) for the Wick-ordered model as a function 
of the separation between the planes in dimension D = 3. 

Fig. 3 - Renormalized mass for the non- Wick-ordered model as a function of the spacing L 
between the planes 

Fig.4 - Renormahzed couphng constant (in units of j^) for the non- Wick ordered model as a 
function of the distance between the planes in dimension D = 3. 
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